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^ \ Abstract 
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^ ' The Ising model in the presence of a random field is investigated within the mean 

field approximation based on Landau expansion. The random field is drawn from 

C^ \ the trimodal probability distribution P{hi) = p6{hi — Hq) + q6{hi + Hq) + r5{hi), 

c/3 ' where the probabilities p, q, r take on values within the interval [0, 1] consistent 

with the constraint p + q + r = 1 (asymmetric distribution), hi is the random 

field variable and /iq the respective strength. This probability distribution is an 

extension of the bimodal one allowing for the existence in the lattice of non magnetic 

particles or vacant sites. The current random field Ising system displays second order 



Q . phase transitions, which, for some values of p,q and /iq, are followed by first order 

CJ I phase transitions, thus confirming the existence of a tricritical point and in some 

cases two tricritical points. Also, reentrance can be seen for appropriate ranges 
of the aforementioned variables. Using the variational principle, we determine the 
equilibrium equation for magnetization, solve it for both transitions and at the 
lO ' tricritical point in order to determine the magnetization profile with respect to h^. 
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1 Introduction 



The pure models for crystalline materials can describe the respective exper- 
imental samples in exceptional situations, since such a sample can contain 
impurities, broken bonds, defects, etc., making real physical systems to never 
be translationally invariant, thus necessitating the modification of pure models 
appropriately for comparing the experimental results with the theoretical pre- 
dictions. A small amount of quenched randomness can infiuence significantly 
the phase transitions replacing a first-order phase transition (FOPT) by a 
second-order phase transition (SOPT), so that tricritical points and critical- 
end-points are suppressed [1]. In two dimensions, an infinitesimal amount of 
field randomness destroys any FOPT. One such situation is the presence of 
random magnetic fields acting on each spin in an otherwise free of defects 
lattice; the respective pure system is considered to be described by an Ising 
model, which is now transformed into the random field Ising model (RFIM) 
in the presence of random fields [2,3,4]. RFIM had been the standard vehicle 
for studying the effects of quenched randomness on phase diagrams and crit- 
ical properties of lattice spin-systems and had been studied for many years 
since the seminal work of Imry and Ma [4]. Associated with this model are 
the notions of lower critical dimension, tricritical points, higher order critical 
points and random-field probability distribution function (PDF). The simplest 
model exhibiting a tricritical phase diagram in the absence of randomness is 
the Blume-Capel model - a regular Ising spin-1 model [5,6]. Although much 
effort has been invested for the study of the RFIM, the only well-established 
conclusion is the existence of a phase transition for rf > 3 (d space dimen- 
sion), that is, the critical lower dimension di is 2 after a long controversial 
discussion [4,7], while many other issues are still unanswered; among them is 
the order of the phase transition, the existence of a tricritical point (TCP) 
and the dependence of these on the form of the random field PDF. Accord- 
ing to the mean-field approximation (MFA) the choice of the random field 
PDF can lead to a continuous ferromagnetic/paramagnetic (FM/PM) bound- 
ary as in the single Gaussian, whereas for the bimodal this boundary can 
be divided into two parts, an SOPT branch for high temperatures and an 
FOPT branch for low temperatures separated by a TCP at kTl/{zJ) = 2/3 
and hi/{zJ) = {kT^^/{zJ)) x argtanh(l/v/3) ~ 0.439 [8,9,10], where z is the 
coordination number, k is the Boltzmann constant and T^*, /i* are the tricrit- 
ical temperature and random field, respectively, such that for T < T^ and 
h > hi the transition to the FM phase is of first order. However, this behav- 
ior is not fully elucidated since in the case of the three dimensional RFIM, 
the high temperature series expansions yield only continuous transitions for 
both PDF's [11]; according to Houghton et al [12] both distributions predict 
a tricritical point with /i* = 0.28 ± 0.01 and Tj: = 0.49 ± 0.03 for the bimodal 
and al = 0.36 ± 0.01 and T^ = 0.36 ± 0.04 for the Gaussian with critical 
standard deviation a*. Galam and Birman studied the crucial issue for the 



existence of a TCP within the mean-field theory for a general PDF p{H) by 
using an even-degree free energy expansion up to eighth degree in the order pa- 
rameter; they proposed some inequalities between the derivatives of the PDF 
up to sixth order at zero magnetic field for the possible existence of a TCP 
[13]. In Monte Carlo studies for d = 3, Machta et al [14], using the Gaussian 
distribution, could not reach a definite conclusion concerning the nature of 
the transition, since for some realizations of randomness the magnetization 
histogram was two-peaked (implying an SOPT) whereas for other ones was 
three-peaked implying an FOPT; Middleton and Fisher [15], using a similar 
distribution for T = 0, suggested an SOPT with a small order-parameter 
exponent /3 = 0.017(5); Fytas et al [16], following Wang-Landau and Lee en- 
tropic sampling schemes for the bimodal distribution function with ho = 2 and 
ho = 2.25 for a simple cubic lattice, concluded that their results indicated an 
SOPT by applying the Lee-Kosterlitz free energy barrier method; Hernandez 
and coworkers claim they have found a crossover between an SOPT and an 
FOPT at a finite temperature and magnetic field for the bimodal distribution 
function [17]. One of the main issues was the experimental realization of ran- 
dom fields. Fishman and Aharony [18] showed that the randomly quenched 
exchange interactions Ising antiferromagnet in a uniform field H is equiva- 
lent to a ferromagnet in a random field with the strength of the random field 
linearly proportional to the induced magnetization. Also another interesting 
result found by Galam [19] via MFA was that the Ising antiferromagnets in a 
uniform field with either a general random site exchange or site dilution have 
the same multicritical space as the random-field Ising model with bimodal 
PDF. 

The usual PDF for the random field is either the symmetric bimodal 

P{h,) = p6{hi - ho) + q6{hi + ho) (1) 



where p is the fraction of lattice sites having a magnetic field ho, while the 



rest fraction has a field {—ho) and p 
or double symmetric. 



g = I [8,20,21], or the Gaussian, single 
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with mean value zero and {ho,— ho), respectively, and standard deviation cr 

[9,22]. 



Galam and Aharony, in a series of investigations, presented a detailed analysis 



via mean field and renormalization group of a system consisting of n— component 
classical spins (finally choosing n = 3) on a d— dimensional lattice of a uni- 
axially anisotropic ferromagnet in a longitudinal random field extracted from 
a symmetric bimodal PDF {p = q = 1/2) without and with a uniform mag- 
netic field along the easy axis, respectively [23,24]. The uniaxial anisotropy 
was chosen to be along the easy axis and the exchange couplings were of the 
form J*^^) = aJ^^\ where a is the anisotropy and < a < 1. Depending on 
the anisotropy (small, medium, large) a variety of phases (longitudinal, trans- 
verse, paramagnetic), critical, bicritical, critical-end points as well as a mul- 
ticritical point (an intersection of bicritical, tricritical and critical-end-point 
lines) resulted. In addition to these purely theoretical investigations, Galam, 
proposed a model (diluted random field) in his attempt to reproduce some 
of the features in the phase diagram of the experimental sample consisting of 
the mixed cyanide crystals X{CN)xYi_,ji., where X stands for an alkali metal 
(K,Na,Rb) and Y a spherical halogen ion (Br, CI, I); the dilution of the pure 
crystal XCN is achieved by replacing CN by the halogen ions Y [25]. The 
pure alkali-cyanide XCN crystal ferroelastic transition disappears at some 
concentration x^ of the cyanide; its numerical value depends on both compo- 
nents X, Y . By choosing a model Hamiltonian (ferromagnetic Ising-type with 
nearest neighbor interaction) with dilution and a symmetric trimodal PDF 
for the random fields, Galam, using MFA, managed to predict the involved 
first and second order phase transitions with the interfering TCP as well as 
the respective concentration for a phase transition to occur depending on the 
procedure considered. The random fields were necessary because there were 
experimental evidences that below x^ cyanide displayed orientational freezing 
and the random fields were used for fixing this orientation. The involved prob- 
ability pt in PDF as well as the critical threshold Xc were expressed in terms 
of microscopic quantities. 

Recently, the asymmetric bimodal PDF (1) with p 7^ g, in general, has also 
been studied in detail [26]. This study has revealed that for some values of 
p and /iq the PM/FM boundary is exclusively of second order; however, for 
some other ranges of these variables this boundary consists of two branches, 
a second order one and another of first order, thus confirming the existence 
of a tricritical point, whose temperature depends only on the probability p in 
(1). In addition to these findings, the occurrence of reentrance has been cor- 
roborated as well as complex magnetization profiles with respect the random 
field strength h^. For p = q, symmetric bimodal PDF, the results found by 
Aharony were confirmed [8]. 

An immediate generalization of the asymmetric bimodal (1) is the asymmetric 
trimodal one, 

P{h,) = p6{h, - ho) + q6{h, + ho) + r6{h,) (3) 



where p+q+r = 1. In earlier studies, the probabihties p, q had been considered 
as equal and related to r by the relation p = q = {^ — ''")/2, symmetric 
PDF [27,28]. The third-peak, introduced in addition the other two ones in 
the bimodal (1) and associated with the third term in (3), is to allow for the 
presence of non magnetic spins or vacancies in the lattice that are not affected 
by the random magnetic fields and reduces the randomness of the system, as 
well. 

For the critical exponents of the three-dimensional RFIM, it seems that there 
is broad consensus concerning their values except for the specific heat exponent 
a, for which there is much dispute concerning its numerical value, since its 
sign is widely accepted to be negative. The main source of information for the 
critical exponents are Monte Carlo simulations. However, they provide various 
values depending on the probability distribution considered. Middleton and 
Fisher concluded that the a-exponent is near zero, a = —0.01 ± 0.09 [15]. 
Rieger and Young, considering the bimodal distribution, found a = — 1.0 ± 
0.3 [29], Rieger, using the single Gaussian distribution, found a = —0.5 ± 
0.2 [30], whereas Hartmann and Young, in ground-state calculations, found 
a = -0.63 ± 0.07 [31]. Nowak et al concluded that a = -0.5 ± 0.2 [32], 
whereas Dukovski and Machta found a positive value, namely, a = 0.12 [33]. 
Malakis and Fytas [34], by applying the critical minimum-energy subspace 
scheme in conjunction with the Wang-Landau and broad-histogram methods 
for cubic lattices, proved that the specific heat and susceptibility are non- 
self-averaging for d = 3 using the bimodal distribution. The same ambiguous 
situation prevails in experimental measurements, see Ref. [35]. 

In this work, we study the RFIM with the asymmetric trimodal PDF (3) with 
arbitrary values for the probabilities p, q {p + q + r = 1) in order to investigate 
the phase diagrams, phase transitions, tricritical points and magnetization 
profiles with respect to ho. The paper is organized as follows: In the next 
section, the suitable Hamiltonian is introduced and the respective free energy 
and equation of state for the magnetization are derived. In section 3, the phase 
diagram, tricritical points and magnetization profiles for various values of p 
and q are calculated and discussed; we close with the conclusions in section 4. 

2 The model 

The pure Ising model Hamiltonian in the presence of random fields changes 
into, 

H = -jJ2S,S,-J2h,S,, S, = ±l (4) 

The summation in the first term extends over all nearest neighbors and is 
denoted by < i,j >; in the second term hi represents the random field that 



couples to the one-dimensional spin variable Si. We also consider that J > 
so that the ground state is ferromagnetic in the absence of random fields. The 
presence of randomness necessitates considering two averaging procedures, 
the usual thermal average, denoted by angular brackets (...), and the disorder 
average over the random fields denoted by {■■■)h for the respective PDF. We 
also make assumptions concerning the random field hi, 

< hi >h= {p - q)ho, < hihj >h= hl6ij (5) 



The former relation in (5) vanishes for a symmetric PDF {p = q), whereas 
for the asymmetric PDF {p ^ q) it is non zero implying that the system has 
residual magnetization, thereby affecting considerably the system's magneti- 
zation; a similar case has appeared in Ref. [26]. The latter one implies that 
there is no correlation between hi at different lattice sites. 

According to the MFA, the Hamiltonian (4) takes the form [8,9,20,26], 

Hmfa = ^NzJM^ - J2{zJM + h,)S, (6) 



where N is the number of spins in the lattice and M the magnetization; the 
respective free energy per spin within the MFA is. 



j^{F)h = ^zJM^ - ^ (ln{2cosh[/3(z JM + h,)]})H 

= ^zJM^ - ^ /" P{hi) ln{2 cosh[f3{zJM + hi)]}dhi (7) 



where the probability P{hi) is chosen to be the trimodal (3) and p, q take on 
any value within the interval [0,1], consistent with the relation p + q + r = 1, 
and (5 = l/{kT). 

The magnetization is the solution to the equation d{{F)h/N)/dM = 0, equi- 
librium condition, 

M = {tanh[/3{zJM + hi)])h (8) 

If the distribution P{hi) is symmetric, P{hi) = P{—hi), which occurs for 
p = q = {1 — r)/2, then the value M = (PM phase) will always be a solution 
to (8), otherwise this is not if P{hi) is non symmetric, p ^ q] however, this 
can be remedied if an auxiliary field Vq is introduced into the system such that 



(tanh[/3(/i, + Vo)])h = (9) 



inducing, in this way, the PM phase. However, this relation acts also as a 
constraint on the system under consideration influencing, nevertheless, its be- 
havior. The free energy (7), in the presence of the auxiliary field Vq, takes now 
the form. 



l^{F)h = yjM^-^{\n{2cosh[/3{zJM + h^ + Vo)]})H 
2 /3 •- ° 2! 3! 4! 



6! 



M^} (10) 



after expanding the quantity in angular brackets in powers of M and calculat- 
ing the average values using (3) with a = f3Jz. By setting ti = tanh[/3(Vo+/ii)], 
t+ = tanh[/3(Vo + ho)], t_ = tanh[/3(Vo — ^0)] and to = tanh[/9V(3], we have, 

Fo=(ln{2cosh[/3(K) + M})/^ 

= ln2+plncosh[/3(\/o + ho)] + gin cosh [/3(\/o - ho)] + r lncosh[/3\/o] 
Fi = {ti)h = pt+ + gt- + rto 
F2 = {1- tl)h = l-ptl- qtl - rtl 
Fs = {-2U{l-tf))h 

= -2pt+(l - tl) - 2gt_(l - ti) - 2rto(l - 4) 
F, = {2{l-tf){3t^-l)), 

= 2p{l - tl){3tl - 1) + 2g(l - tl){3tl - 1) + 2r(l - tl){3tl - 1) 
F, = {8il-t^)il5tt-15t^ + 2)), 

= 8p{l - tl){l5tl - 15tl + 2) + 8g(l - t^)(15ti - 15t^ + 2) (11) 

+8r{l-tl){154-15tl + 2) 

The condition (9) for the existence of the PM phase is equivalent to setting 

pt+ + qt^ + rto = (12) 

The equilibrium condition d{{F)h/N)/dM = yields, 

M = aF,M + ^M' + ^M' + ^M' (13) 

2! 3! 5! ^ ' 



or, 



M = AM + BM^ + CM^ + EM^ (14) 



A = aFo,B = 



2! 



,C^ 



a^F^ ^ _ a^Fa 



3! 



.E = 



5! 



(15) 



In RFIM if there is a phase transition it will be associated with the mag- 
netization and the possible two phases are the PM with M = and FM 
with M 7^ 0. The phase boundary is found by solving Eq. (14) in conjunc- 
tion with the free energy (10) in case of an FOPT. The SOPT boundary is 
determined by setting A = 1 and C < 0, whereas the FOPT boundary by 
A = 1 and C > 0. These two boundaries, whenever appear sequentially, are 
joined at a tricritical point determined by the condition A = 1 and C = 
[8,20,21,22,26,36,37,38], provided that E < (equivalently, Fq < 0) for sta- 
bility, as also in [36,37]. However, for the FOPT boundary we shall use, in 
addition to (12) and (14), the requirement of the equality of the respective 
free energies, F{M = 0) = F(M ^ 0), where F = {F)h/N. 

3 Phase diagram. Magnetization profiles 
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Fig. 1. The tricritical temperature against the probabihty q for specific values of 
p. In panel (a) {p = 0.01) it varies non-monotonically with a maximum value at 
q = 0.50. In panel (b) (p = 0.40) the left-hand-side points correspond to tricriti- 
cal temperatures for low q- values (0.0 < q < 0.15). The right-hand-side group of 
points forms two branches, the upper one refers to the upper tricritical tempera- 
tures (0.35 < q < 0.60), whereas the lower one corresponds to the lower tricritical 
temperatures for fewer values of q, 0.35 < q < 0.46. In panel (c) (p = 0.48) the 
tricritical temperature varies non-monotonically with a minimum at g = 0.25. 

Using the conditions for the calculation of the TCP together with Eq. (12), the 
TCP coordinates (T^*"^, h^'"^, V^'^^) are calculated as functions of the prob- 
abilities p, q. These exist only for a limited number of p's and g's, namely, p G 
[0.0, 0.63] whereas the respective g— values depend on the specific p— values, 
but they lie in same interval, as well. However, the tricritical temperature T^^^ 
does not satisfy any simple closed-formula as in the asymmetric bimodal PDF 
[26]. The resulting tricritical temperatures exhibit a variety of variations as 
functions oi p,q, see Fig. 1. However, for some p- and q-values two tricritical 
temperatures (upper and lower ones) occur, see Fig. 1(b) [39,40], since both 
temperatures are solutions to the simultaneous equations aF2 = 1, F4 = 0. 
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Fig. 2. Tricritical point coordinates for the bimodal PDF, (a) tricritical temperature 
and (b) random field /iq and auxiliary field Vq, resulting form trimodal PDF in the 
limiting case r = 0. The agreement of these figures with the corresponding ones in 
Ref. [26] is complete. 
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Fig. 3. Variation of the random field strength /iq and auxiliary field Vq with q for 
specific p— values at the tricritical point for the trimodal PDF. The random field 
/iQ exhibits a non-monotonic behavior in (a) for p = 0.01 and monotonic in (b) for 
p = 0.50. The auxiliary potential Vq increases monotonically with q for any p. Both 
quantities are in units of (Jz), i.e., /iq = Iiq/^Jz), Vq = Vq/{Jz). 



In order to examine the validity of the process under consideration, we focus 
on the asymmetric bimodal PDF resulting from the trimodal by setting r = 
and studied earlier [26]; using the data for r = 0, we recover exactly the same 
plots for the tricritical temperature as well as ho, Vq, see Fig. 2. 

Now, using the relations A = 1 and C = (equivalently, aF2 = 1 and F4 = 0) 
together with (12), the remaining coordinates Hq^'^ and V^*"^ for the TCP 
are calculated, see Fig. 3. Both quantities display two modes of variation; in 
mode-1 {p < 0.49) /i^*"^ varies non-monotonically with q, Fig. 3(a); in mode-2 
{p > 0.50) Hq'"^ decreases monotonically with q, Fig. 3(b). However, in both 
modes, V^^^ increases monotonically, but for type-1 the increase is more steep 
than in type- 2. 

The magnetization at the TCP is found by solving Eq. (13) taking into con- 
sideration the appropriate conditions for the TCP, 

^M^ + ^M^ = (16) 

2! 5! ^ ' 



or 

Fqco^ + eOF^co^ = (17) 

where u = aM, whose solutions are, 

a;r^ = (18) 

u;r^ = (-60Fs/F,f" (19) 



from which the magnetizations M'[2^ = uJi2^*{kT'^'"^/{Jz)) can be deduced. 
The first one, (18), is the magnetization of the PM phase [Mf'-''^ = 0), whereas 
the second one, (19), is the magnetization of the FM phase (Mj'^'^ 7^ 0). 
However, the nonzero solution Mj'-^'^, in general, has lower free energy than 
the zero solution (18), implying that this is the stable solution at the tricritical 
point, see Fig. 4. According to the first relation in Eqs. (5), for the general case 
p ^ q, the mean value of the random field is non zero; this is equivalent to the 
presence of an external magnetic field in the system so that the magnetization 
at the tricritical point scales as Mt = M{T = T'^'^^) ~ hJcp, where h^cp is 
the magnetic field and the tricritical exponent St = 5 according to the Landau 
theory [41,42,43]. In case of equal partial probabilities {p = q), the M^"'^-^- 
magnetization vanishes (Mj*-^-^ = = M'f'~'^) and the system has only the 
zero solution (double root), which now becomes the stable one. Representative 
plots of the magnetization M'^^^ with respect to q appear in Fig. 5 for specific 
values of p. 
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Fig. 4. Free energy of the solutions (18) and (19) at the tricritical point for p = 0.55; 
the upper plot (i) corresponds to the zero solution Mi and the lower (ii) to the 
nonzero solution M2. For this p— value, the M2 is the stable solution, whereas the 
Ml is metastable; this happens, in general, for other |)— values. The M2 solution 
coincides with the zero one (Mi) only for p = q. 

In a previous communication [26], the PDF of the RFIM was selected to be 
the asymmetric bimodal (1); this system displayed a symmetric behavior at 
the tricritical point with respect to the probability p ; especially, two distinct 
tricritical points with probabilities pi and p2, respectively, such that pi + 
P2 = 1, have identical tricritical temperatures and random fields, whereas the 
respective auxiliary fields and non zero magnetizations are absolutely equal. 
A similar behavior is also displayed by the present model with respect to the 
probabilities p and q; if the probabilities (^1,^1) and {p2,q2) of the tricritical 
points of two systems are interchanged, namely, P2 = (11,12 = Pi, then these 
systems have the same tricritical temperatures and random fields, whereas 
the respective auxiliary fields V^^^ and the nonzero magnetizations M'^'"^ 
are absolutely equal. These systems also have equal the respective free energies 
for the zero magnetization {F{pi,qi,M'['^'^ = 0) = F{p2,q2,M'['-^^ = 0)) as 
well as the nonzero one {F{pi, gi, Mj*-^-^) = F{p2, q2, ^2^^))-^ the latter result 
implies that the two magnetizations M2'"^{pi,qi), M2^^{p2,q2) are equally 
probable, an expected result, since the magnetizations have equal absolute 
values, the only difference being in their sign, implying that no direction is 
favored. 

After the calculation of the tricritical point coordinates {T'^'-''^, h'^'"^, Vq^*"^) 
for the respective values of p and q, we proceed to determine the phase diagram 
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Fig. 5. Variation of the nonzero magnetization M2 = ( ~ SOFs/Fgj ' with q 
for specific p— values at the tricritical point. M2 is, in general, the stable solution. 
Except for p = 0.01, in the other panels the TCP-magnetization displays a non 
monotonic behavior as a function of q. 

by solving Eq. (14), taking into consideration the respective conditions for the 
FOPT and the SOFT. By varying the parameters p and q many different 
types of phase diagrams result, as those appearing in Fig. 6(a) labelled by 
the individual p and q values. The curves in Fig. 6(a) are classified into two 
main groups: the first one includes those curves that do not possess a TCP 
corresponding to an SOFT only, curves (i) and (ii); the second group includes 
those having a TCF, which joins the FOFT-branch with the SOFT-branch 
of the phase diagram, curves (iii) and (iv). The occurrence of an FOFT, and 
subsequently of a TCF, results from the competition between the first term in 
the Hamiltonian (4) (tending to make parallel the spins) and the second term 
of random forces inducing disorder. For small values of h^ the competition is 
weak allowing the first term to dominate, but as h^ increases the second term 
dominates over the first one, thus changing the phase transition from second 
order to first order. Some of the curves in the latter group have a second 
tricritical point and/or present reentrance Fig. 6(a(iv)). Reentrance might also 
be attributed to the competition between these two terms in the Hamiltonian 
(4). In the phenomenon of reentrance a vertical line in the (/iq, T)-plane crosses 
the transition line twice, in that, by lowering the temperature at constant /iq, 
one observes first a PM/FM transition and then, on further lowering the 
temperature, an FM/PM transition appears so that the magnetization is zero 
although the temperature is low and the system remains in the FM phase for 
these temperatures. Fig. 6(b(i)) or another transition may take place from 
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Fig. 6. Phase diagram of the Hamiltonian (4) in the (ho — T) plane for specific values 
of p and q. The solid curve is a line of critical points and the dashed one is a line of 
first-order phase transitions, joined smoothly by a tricritical point (full circle). The 
system has, in general, second-order transition as in panel (a) for p = 0.35, q = 0.0, 
curve (i) and p = 0.35, q = 0.20, curve (ii); however, for p = 0.05, g = 0.40, curve 
(iii), it displays both phase transitions, FOPT and SOPT joined by a tricritical 
point, whereas for p = 0.45, g' = 0.45, curve (iv), the system has two tricritical 
points, occurring also for other values of p and q. In the latter case, the system 
displays reentrance, as well. In panel (b) the phenomenon of reentrance is shown in 
enlargement; for p = 0.45, q = 0.50 (curve (i)) the system remains in the PM phase 
for low temperatures and medium random fields, whereas for p = 0.45, q = 0.55 
(curve (ii)) it returns to the FM phase for low temperatures and higher random 
fields. Temperature T is expressed in units of (Jz/k), i.e., T = kT/{Jz). 

PM to FM with the system, now, in the FM phase for low temperatures 
and high fields, Fig. 6(b(ii)); occasionally, the region of the FM phase shrinks 
significantly, Fig. 6(a(iv)). The vanishing of magnetization for high values of 
p (j9 ~ 0.45, 45% up-spins) and small q values (where one would expect a 
nonzero magnetization) can also be due to the presence of the auxiliary field 
Vo, which annihilates the excess magnetization making the system to behave 
like an antiferromagnet. This phenomenon is evident from the bending of the 
phase transition lines on lowering the temperature thus forming an inverted 
"C" ("boomerang" shape) and appearing in enlargement in Fig. 6(b). The 
respective p values lie in the interval [0, 0.50], whereas the q values lie in a much 
smaller interval, namely, [0.40, 0.55]; the only exception is for p = 0.45 when q 
takes on values in the interval [0,0.55]. However, within the MFA reentrance 
may lead to nonphysical values (negative) for the specific heat, since energy 
will also present reentrant behavior as magnetization because energy, within 
MFA, is proportional to the magnetization squared thus behaving similarly. 
The two successive transitions can be of first or second order depending on 
p, q and ho. 
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Solving Eq. (14) to determine the phase diagram, the magnetization is also 
calculated for either phase transition. The condition A = 1 or aF2{f3, Vq, h^) = 
1 leads to 

ptl + qtl+rtl = ^^^ (20) 

and, by setting T2 = pt"^ + qt^__ + rtp, Eq. (20) can be written as, 

T2 = ^^ (21) 

a 

Inverting Eq. (21) the respective temperature for either phase transition can 
be determined, 

f^l-T, (22) 

In order to specify the type of the transition, the sign of C = 0:^^4/6 is 
checked; however, to facilitate the calculations, the quantity C is rewritten as, 

C = ^[4T2 - 3T4 - 1] = tf^[l - T4 - ;i] (23) 

using Eq. (21) and T4 = pt^ + gti +r tp. For an SOFT, C is negative [8,9,26], 
then Eq. (23) yields, 

Ta>1-^ (24) 

6a 

otherwise if 

T4 < 1 - — (25) 

6a 

the resulting transition is an FOFT. However, in order to determine the mag- 
netization for an FOFT the expression (13) is combined with the equality of 
the respective free energies, 

F(M = 0) = F(M ^ 0) (26) 

or, 

M2 = F^aM^ + ^a^M^ + ^a^'M* + ^a^M^ (27) 

3 12 360 ^ ' 

Combining Eqs. (13), (27) and using the condition aF2 = 1, leads to, 

Fqco^ + IOF4W = (28) 



which is split up into two equations, the first is cji = or, equivalently. Mi = 0, 
FM phase, whereas the other is, 

Fecu^ + IOF4 = (29) 
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leading to the nonzero solutions, FM phase, 



UJ2,3 = ±V-10F4/F6 



(30) 



for which Fq < for stability requirements as far as F4 > for an FOPT; the 
value for F3, consistent with (13) and (27), is F3 = — (^4/6)1/— IOF4/F6 for the 



positive root in (30) and F3 = (^4/6)1/— IOF4/F6 for the respective negative 
root. From the solution of this equation we can extract the magnetizations 
M2,3, since the temperature {kT/zJ) is already known from (22). 
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Fig. 7. Magnetization profile vs. /iq, for (a) p = 0.25, g = 0.55, point A is a critical 
point, (b) p = 0.45, q = 0.55 with three critical points A, B, C. Point A is a normal 
critical point. B a critical-end-point; the two critical phases are on either branch 
at point B coexisting with the non critical phase on the branch CD. C a double 
critical-end-point; the first group of the two critical phases are on the branches 
CGA and CE coexisting with the non critical phase on the branch BF, the second 
group of two critical phases are on the branches CD and CB coexisting with the non 
critical phase on the branch BF. This plot represents a closed-loop phase diagram 
with left-hand and right-hand-side critical points. 

Considering, now, the randomness strength h^ as a control parameter, sim- 
ilarly as temperature, we study the variation of the non-zero positive mag- 
netization M2 = U2 * {kT/{Jz)) for an FOPT, Eq. (30), coexisting in equi- 
librium with the zero magnetization Mi = 0, by calculating the respective 
magnetization profile as a function of ho for specific values of p and g; the 
negative M3 = — M2 behaves analogously. These profiles appear in Fig. 7 as 
functions of h^ and look like the ones when the temperature is the control 
parameter for an FOPT. For p = 0.25, q = 0.55 the magnetization has a nor- 
mal critical point Fig. 7(a), whereas for p = 0.45, g = 0.55 this structure is 
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transformed into a closed loop magnetization profile, closed miscibility gap, 
Fig. 7(b); it possesses multiple critical points, the point A is a pure critical 
point, whereas S is a critical-end-point (CEP) because the thermodynamic 
states along the branches BF and BC become identical at the point B in the 
presence of the third non critical phase (spectator phase) along the branch 
CD for the same value of ho; the CEPs can be found in binary fluid mix- 
tures, superfluids, binary alloys, liquid crystals, ferromagnets, ferroelectrics, 
etc. [26,44,45,46,47,48,49,50,51]. A simple example for the occurrence of the 
CEP is to consider a binary fluid mixture with the proper thermodynamic 
conditions so that three fluid phases can coexist in equilibrium and enclosed 
in a capsule; these phases are called V (vapor, top phase), Li (liquid phase 
rich in species 1, middle phase) and L2 (liquid phase rich in species 2, bottom 
phase). The top and middle phases are separated by the interface (VLi), as 
well as the middle and the bottom ones by the interface (L1L2); if the thermo- 
dynamic conditions are such that the interface VLi disappears with the phases 
V, Li becoming identical in the presence of the non critical L2-phase, then 
that thermodynamic state is a CEP at the temperature Tyf_^ and the new 
phase VLi coexists with the L2 one in equilibrium. However, if we consider 
that the interface L1L2 disappears with the phases Li, L2 becoming iden- 
tical in the presence of the non critical V-phase, then that thermodynamic 
state is a CEP at the temperature Tfl^'f^ and the new phase L1L2 coexists 
with the V one in equilibrium. Another occurrence of a CEP is the double 
critical end point (DCEP) where two critical lines end simultaneously at a 
flrst order phase boundary; such a situation appears at the point C, where 
the two critical phases on the branches CCA and CE become identical at 
the point C in the presence of the non critical phase on the branch BF; the 
second group of two critical phases are on the branches CD and CB that be- 
come identical in the presence of the non critical phase on the branch BF; 
these are also observed in binary fluid mixtures, metamagnets, Ising antifer- 
romagnets with next-nearest-neighbor interactions, RFIM, three dimensional 
antiferromagnetic spin-1 Blume-Capel model [26,45,51,52]. Point A can also 
be considered as a left-hand side critical point and 5, C right-hand side ones. 

Now, we consider the values of p and q for which the system exhibits only an 
SOFT; Eq. (14) takes the form for A = 1, 

Fqu^ + 20^4^2 + mF^oj'^ = (31) 



The value ui = (two-fold) is again a solution or, equivalently. Mi = (PM 
phase); the other three roots are the solutions to the third degree equation, 

Fqco^ + 2OF4UJ + 6OF3 = (32) 

which, depending on the values oi p,q and ho, can have either only one real 
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noil zero solution if A = g^ + r^ > (r = -SOFg/Fe, q = {20F4)/{3Fe)), 



UJ2 = \/r + VA + \/r-VA (33) 

or three real non zero ones for A < 0, 

U2 = 2^ cos(^/3) 

00^ = -^ [cos(^/3) + VS sin(^/3)] 

cj^ = -s^ [cos(^/3) - y3sin(^/3)] (34) 



where p = ^/r^ - A, 9 = arctan(y^A/r) and Mi = uJi*{kT/{Jz)), i = 2, 3, 4. 
The solutions for an SOPT are classified into two groups, group-1 includes the 
zero solution (Mi = 0) and the single nonzero one M2 of (33), whereas group- 
2 includes again the zero solution and the nonzero ones M2, M^, M4^ of (34). 
Depending on the value of p, q and ho, there can be transitions between these 
two groups. The majority of the SOPT solutions belong to the group-2; for 
given values of p and q the stable solution is the zero (PM phase) for small 
values of ho, whereas for higher ones the most probable to be stable is the M2 
solution. For the group-1, the zero solution is stable for small values of ho and 
the M2 solution for higher values. 

The investigation was also extended towards the zero-temperature case, T = 0; 
in this case the free energy (7) reduces to. 



^(F)^ = ^zJM^ - i (ln{2cosh[/3(z JM + /i,)]})/. 

= -zJM^ -{\zJM + h\)h 

I 
= -zJM^ - p\zJM + ho\- q\zJM - ho\ - r\zJM\ (35) 



the external potential was omitted. Applying the equilibrium condition dF/dM 
to (35) we get, 

M J'JM + ho\ \zJM-ho\ \zJM\ 

^^^"^ zJM + ho ^^ zJM-ho ^ zJM ^^^^ 



Analyzing Eq. (36) we find a variety of solutions because of the greater number 
of degrees of freedom than in the case of the bimodal PDF [26]. The solution 
M = 1 is a stable one for p + r > ho/zJ, whereas for p + r < ho/ zJ the stable 



1— r ' 
2 ' 

the results found by Sebastianes and Saxena [28] are recovered, that is, the 



one is M = 1 — 2g. If we consider the symmetric trimodal PDF {p = q = ^^) 
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former result (M = 1) is stable for ^-^ > h^/zJ, whereas the latter (M = r) 
for ^^ < ho/zJ, using the current notation. Across the line h^/zJ = p + r in 
the {p, /io)-plane a first-order phase transition occurs between the two ordered 
phases with M = 1 and M = 1 — 2g. In addition to the aforementioned 
two solutions, there are more ones; the result M = 2{p + q) — 1 is stable for 
p — r > ho/zJ, M = 2p — 1 for p ~r < h^/ zJ , M = 1 — 2p for r — p > h^/zJ, 
M = l—2{p+q) for r—p < ho/zJ and M = —1 ior p+r+ho > 0/zJ. In the first 
case, M = 1, the condition p > {ho/{zJ)) implies that the exchange interaction 
J is much stronger than randomness ho, their ratio is always smaller than one, 
thus forcing the system's spins to order according to the first term in (4). The 
alternative condition p < {ho/{zJ)) implies, now, that randomness is not any 
more negligible but strong enough to influence significantly the spins enforcing 
a p-fraction of them to point up and a g-fraction down, to randomly align with 
the local fields, thus, practically, it dominates, so to speak, over the first term 
in Eq. (4) so that M = p — q + r = 1 — 2g. 



4 Conclusions and discussions 

In the current treatment we have determined the phase diagram and discussed 
some critical phenomena of the Ising model under the influence of a trimodal 
random field, an extension of the bimodal one, to allow for the existence of 
non magnetic particles or vacancies in the system, for arbitrary values of the 
probabilities p and q via the mean-field approximation. The competition be- 
tween the ordering effects and the randomness induces a rich phase diagram. 
The system is strongly influenced by the random field, which establishes a new 
competition favoring disorder; this is obvious by the appearance of first-order 
transitions and tricritical points, in addition to the second-order transitions, 
for some values of p and q; the tricritical point temperature has various modes 
of variation as a function of p and q. The trimodal distribution induces reen- 
trant behavior for the appropriate range of p, q and random field ho. For some 
values of p and q the system can be found either in the PM phase or in the 
FM phase for low temperatures and medium and/or high random fields; oc- 
casionally, the part of the phase diagram allocated to the FM phase shrinks 
significantly. A direct consequence of the asymmetric PDF {p ^ q) is the exis- 
tence of residual mean magnetization, a result of < hi >h= {p — (!)ho, making 
the TCP non zero magnetization M2 to be the stable one in comparison to the 
zero one. Mi; however, for p ^ q (symmetric PDF) the residual mean magne- 
tization vanishes as well as the initially TCP non zero magnetization so that 
M2 = Ml, which is now the stable one. Both asymmetric PDFs, bimodal and 
trimodal, confirm the existence of a TCP and, nevertheless, yield similar mag- 
netization profiles as well as reentrance; however, the trimodal one predicts 
also the existence of a second TCP. The tricritical point temperatures for the 
bimodal and trimodal PDF's are independent of the random field strength ho; 
they depend only on the probability p and p, q, respectively. 



Griffiths extended the notion of the critical point to the so-called multicriti- 
cal points, e.g., tricritical, critical-end-point, double critical-end-point, fourth- 
order, ordered critical point, etc. [53]; however, in order to describe these 
points (except the ffist two) the considered expansion of the free energy (10) 
has to be extended to higher-order terms [27,37,54,55] so that the stability cri- 
teria for such a point are satisfied, but this is beyond the scope of the current 
research. 

The Landau theory breaks down close to the critical point (non-classical re- 
gion) because as the transition temperature is approached the fluctuations 
become important and non-classical behavior is observed. A relative criterion, 
called Ginzburg criterion, determines how close to the transition temperature 
the true critical behavior is revealed, or, in other words, it governs the validity 
of the Landau theory close to a critical point [56] . This criterion relies on any 
thermodynamic quantity but the specific heat is usually considered for de- 
termining the critical region around Tc where the mean field solution cannot 
describe correctly the phase transition. The Landau theory is valid for lattice 
dimensionality greater than or equal to the upper critical dimension (i„ = 4 
in case of presence only of thermal fluctuations. However, in the current case 
the presence of random flelds enhances fluctuations causing the critical region 
to be wider than the one due only to the thermal fluctuations [57,58] and the 
upper critical dimension is increased by 2 to du = 6. 

Our results indicate that on increasing the complexity of the model system 
new phenomena can be revealed as in the current case of including asymmetry 
in the PDF; this inclusion induces drastic changes on the phase diagram, such 
as reentrance and two TCPs, thus conflrming the necessity of treating the 
partial probabilities (p, q, r) of the PDF in the most general aspect to get the 
complete phase diagram. A similar situation appears in the model systems in 
Refs. [23,24] wherein the considered complexity has revealed a rich variety of 
phase diagrams with known and new multicritical points. The results obtained 
in the current investigation by using the MFA need further analysis; these can 
provide a basis for a comprehensive analysis by more sophisticated methods as 
well as experimental implementation. However, they are of no less importance, 
since they show, nevertheless, the expected phenomena to be observed. 
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